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ON ■ Abstract 

(— I , We discuss the Meissner effect for a color superconductor 

0^1 formed by cold dense quark matter. Though color and ordinary 

o ' electromagnetism are broken in a color superconductor, there is a 

1^ . linear combination of the photon and a gluon that remains massless. 

Consequently, a color superconducting region may be penetrated by 

an external magnetic field. We show that at most a small fraction 

^ ■ of the magnetic field is expelled, and if the screening distance is the 

smallest length scale in the problem there is no expulsion at all. We 
calculate the behavior of the magnetic field for a spherical geometry 
relevant for compact stars. If a neutron star contains a quark matter 
core, this core is a color superconductor. Our results demonstrate that 
such cores admit magnetic fields without restricting them to quantized 
flux tubes. Such magnetic fields within color superconducting neutron 
star cores are stable on time scales longer than the age of the universe, 
even if the spin period of the neutron star is changing. 



1 Introduction 

Conventional superconductors result from a condensate of Cooper pairs 
of electrons. Because the Cooper pairs have nonzero electric charge, 
electromagnetic gauge invariance is spontaneously broken: the photon gets 
a mass and weak magnetic fields are expelled by the Meissner effect. In 
this paper we show that the same is not in general true for the color 
superconducting state that is formed by cold dense quarks |l|, ||, |^, H, ||, |, , 
even though here also the Cooper pairs have nonzero electric charge. The 
reason is that a color superconductor is not quite an electric superconductor: 
it makes the gluons massive (there is a color Meissner effect) but does 
not simply make the photon massive. Rather, one linear combination of 
the photon and a gluon becomes massive, but the orthogonal combination 
remains massless. Thus a region of color superconductor can allow itself to be 
penetrated by the component of an external magnetic field that corresponds 
to the unbroken generator. As we will see below, in the limit in which the 
screening length is the shortest length scale in the problem, the magnetic 
field within the color superconductor has the same strength as the applied 
external magnetic field. There is no Meissner effect. Though the interior 
field is not diminished in strength, it is "rotated" relative to the external 
field: it is associated with the U{1) symmetry which is unbroken within 
the superconductor, not with the U{1) of ordinary electromagnetism. If 
the penetration length is not smaller than the thickness of the boundary of 
the color superconducting region (to be defined below), there is a partial 
Meissner effect. In this case, the strength of the field which penetrates the 
superconductor depends on details of the geometry, the relative sizes of the 
screening length and the boundary thickness, and the relative strengths of 
electromagnetism and the color force. 

The most likely place to find superconducting quark phases in nature 
is in the core of neutron stars. If neutron stars achieve sufficient central 
densities that they feature quark matter cores, these cores must be color 
superconductors: because there are attractive interactions between pairs 
of quarks which are antisymmetric in color, the quark Fermi surfaces in 
cold dense quark matter are unstable to the formation of a condensate of 
quark Cooper pairs. Present theoretical methods are not sufficiently accurate 
to determine the density above which a quark matter description becomes 
appropriate, and thus cannot answer the question of whether quark matter, 
and hence color superconductivity, occurs in the cores of neutron stars. What 



theory can do is analyze the physical properties of dense quark matter and, 
eventually, make predictions for neutron star phenomenology, thus allowing 
astrophysical observation to settle the question. 

There are a number of avenues which may allow observations of 
neutron stars to answer questions about the presence or absence of color 
superconductivity. (Examples we do not pursue in this paper include analysis 
of cooling by neutrino emission ||^ and analysis of r-mode oscillations [§].) 
Since neutron stars have high magnetic fields (10^ to lO^'^'^ Gauss in typical 
pulsars [0; perhaps as high as 10^^ Gauss in magnetars |jll|) one prerequisite 



to making contact with neutron star phenomenology is to ask how the 
presence of a superconducting core would affect the magnetic field. We 
will see in this paper that the strength of the magnetic field within the 
superconducting core is hardly reduced, and the magnetic flux within the 
color superconductor is not restricted to quantized flux tubes. The latter 
constitutes a qualitative difference between conventional neutron stars and 
those with quark matter cores. We will see in Section 5 that, unlike in 
conventional neutron stars, the magnetic field within a color superconducting 
core does not vary even if the spin period of the neutron star is changing. 

1.1 A fiducial example 

To give the reader some sense for typical scales in the problem, we now 
describe a fiducial example, which we will use in particular in Sect. ^. The 
numbers in this paragraph make the crude assumption that the quarks are 
noninteracting fermions, and so should certainly not be construed as precise. 
Consider quark matter with quark chemical potential /i = 400 MeV, made of 
massless up and down quarks and strange quarks with mass Mg = 200 MeV. 
(Ms is a density dependent effective mass; this adds to the uncertainty in its 
value.) If the strange quark were massless, quark matter consisting of equal 
parts u, d and s would be electrically neutral. In our fiducial example, on 
the other hand, electric neutrality requires a nonzero density of electrons, 
with chemical potential /Ze = 24 MeV. Charge neutrality combined with the 
requirement that the weak interactions are in equilibrium determine all the 



chemical potentials and Fermi momenta: 



fJ'U 


= /x - |/Xe = 384 MeV, 


Pf ~ f^uj 


f^d 


= /i + |/ie = 408 MeV, 
= yu + |/ie = 408 MeV, 


Pf = f^d, 

(1.1) 


l^s 


p^ = ^^l - M2 = 356 MeV, 




/ie = 24 MeV, 


Pf = Pe ■ 


bar 


yon number density ps = 


-- (l/37r2)[(p^)3 + (pfY + (p^)3] is 5 times 



nuclear matter density. 

A variety of estimates suggest that the gaps at the Fermi surfaces resulting 
from quark-quark pairing are about A ~ 20 — 100 MeV, resulting in critical 
temperatures above which the color superconductivity is destroyed which are 
about Tc ~ 10 — 50 MeV. Neutron stars have temperatures T <^ Tc, and if 
they have quark matter cores these cores are certainly in the superconducting 
phase. Similar estimates for A and Tc are arrived at either by using 
phenomenological models, with parameters normalized to give reasonable 
vacuum physics |P, |, ||, ^ |12|, |l^, |l^, |^ or by using weak coupling methods 



|T|, [13, 0, |T|, |20|, gl], g, valid for ^ ^ cx) where the QCD coupling 
g{fi) does become weak. Neither strategy can be trusted quantitatively for 
/i ~ 400 MeV, where g{p) ~ 3, but it is pleasing that both strategies agree 
qualitatively. 

1.2 The CFL and 2SC phases of quark matter 

There are, in fact, two different superconducting phases possible, which have 
very different symmetry properties. If A is large compared to the differences 
among the three Fermi momenta in ( |1.1|) , (ud), (us) and (ds) condensates 
all form. Chiral symmetry is broken by color-flavor locking [^ in this phase. 
This is the favored phase at large /i where the differences between Fermi 
momenta decrease ^^ ^. This CFL phase has the same symmetries as 
baryonic matter which is itself sufficiently dense that the hyperon and nucleon 
densities are all similar, and there need not be a phase boundary between 



CFL matter and baryonic matter |26, 24, 25|. The properties of the CFL 



phase have been further investigated in Refs. [ P^ , 

Now, imagine starting with CFL matter at very large n and reducing 
fj,. Because of the nonvanishing strange quark mass, as /x decreases the 



differences among the Fermi momenta increase. It may happen that before n 
has decreased so far that a quark matter description ceases to be vahd, one 
may find a phase of quark matter in which only two flavors {u and d) and two 
colors (chosen spontaneously) of quarks pair. Chiral symmetry is restored in 
this two-flavor superconductivity (2SC) phase, which is also the phase which 
arises in QCD with no strange quarks at all |l|, 0, 0, ^, ||. Because nature 
chooses a strange quark which can neither be treated as very light nor as 
very heavy, present theoretical analyses are not precise enough to determine 
whether quark matter at densities typical of neutron star interiors is in the 
CFL phase, or whether in this range of n the 2SC phase is favored. For 
example, in (|1 . 1|) the splitting between d and s Fermi momenta is ~ 50 MeV, 
of the order of typical gaps. Current theoretical methods are not reliable 
enough to determine whether the quark-quark interactions in QCD are strong 
enough to generate a (ds) condensate larger than 50 MeV, as required in the 
CFL phase, or are somewhat weaker, admitting only the 2SC pairing. 

In discussing conventional superconductors in a magnetic field, 
one normally begins by distinguishing between Type I and Type II 
superconductors. Color superconductors are Type I at asymptotically high 
densities |1^, ^, Q, but they may be Type I or Type II at neutron 



star densities. However, this distinction will not be of importance. First 
of all, the thermodynamic critical field H^ required to destroy the color 
superconductivity is on the order of 10^^ Gauss. If color superconductors 
exhibited the conventional Meissner effect, whether or not they were Type I 
or Type II they would simply exclude neutron star magnetic fields, which are 
in fact weak for our purposes. Second of all, we shall see that the presence 
of an unbroken rotated electromagnetism changes the story completely.^ 
In the next three sections, we solve the problem of how a sphere of color 
superconducting matter, which can admit a rotated magnetic field, responds 
to an applied ordinary magnetic field. In the final section, we discuss the 
consequences of our findings for neutron stars. 

2 Rotated electromagnetism 



^ See g, |29| for early attempts to analyze magnetic fields in color superconductors. 
These authors neglect the existence of the unbroken rotated electromagnetism. 



2.1 The new photon 

The fundamental fields in the theory are the quarks ipf (with color index a 
and flavor index i) and the gauge flelds: the photon A^ and the gluons G" 
n = 1 . . .8. At low temperature and high density, the quarks form Cooper 
pairs, associated with a Higgs fleld 0^ ~ '4^t4'j- This fleld gets a vacuum 
expectation value which takes the form 

«) ~ Ae-^e,, (2.1) 

in the 2SC phase, with i and j running over u and d only and a and (3 
running over two of the colors only. In the CFL phase, quarks of all three 
flavors and colors pair and the expectation value of the fleld takes the form 
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> ~ AiSfS^ + A^SfSf , (2.2) 



with all indices now taking on three values. In both the CFL and 2SC phases, 
the condensate leaves an unbroken f/(l) generated by 

Q = Q + r]Ts, Q(0) = O, (2.3) 

where rj = l/y/S for CFL, and r] = — 1/(2-\/3) for 2SC. Q is the conventional 
electromagnetic charge generator, and Tg is associated with one of the gluons. 
In the representation of the quarks, 

Q = diagd, — |, — |) in flavor u, d, s space 

(2.4) 
Tg = 4s diag(l, 1, —2) in color r, g, b space. 

As is conventional, we have taken tr(T8T8) = 2. The Q"Charge of all the 
Cooper pairs which form the condensates are zero. 

To see exactly which gauge fleld remains unbroken, look at the covariant 
derivative of the Higgs fleld: 

D,{<P) = [d, + eA^QW + gGlT^'^^) (0) (2.5) 

From ( |2.3| ), we see that the kinetic term |i5(0)p will give a mass to one gauge 
fleld 

^jri^e^ + g^ 



but the orthogonal hnear combination 

AQ = ''' I = cos aoA, + sin aoGj (2.7) 

y T^^e^ _j_ gZ 

will remain massless. The denominators arise from keeping the gauge field 
kinetic terms correctly normalized, and we have defined the angle ao, 

cosao = —= , (2.8) 

^^ZgZ _|_ gZ 

which specifies the unbroken f/(l). At neutron star densities the gluons are 
strongly coupled (g'^/^An) ~ 1), and of course the photons are weakly coupled 
(e^/(47r) ^ 1/137), so cto — Tje/g is small: ao ~ 1/20 in the 2SC phase and 
ao ~ 1/40 in the CFL phase. In both phases, the "rotated photon" consists 
mostly of the usual photon, with only a small admixture of the Tg gluon. 
We can now see that the electron couples to A^ with charge 



egj yrffp' + g"^ which is less than e, so the new photon is slightly more weakly 
coupled to electrons than the old one. It can also be shown that in the CFL 
phase, the quarks have Q-charges which are integer multiples of the Q-charge 
of the electron . 

2.2 The Meissner effect 

The main purpose of this paper is to study the Meissner effect in a high 
density region, which we will assume to be a sphere, like the core of a neutron 
star, with radius R. In this region, color superconductivity occurs. The 
unbroken generator Q is associated with a rotated electromagnetic field A^, 
with rotation angle ao, as described above. The broken generator X gives no 
long-range gauge field. For simplicity we will treat the region outside the core 
not as nuclear matter, but as vacuum. As we will explain in Sect. |^, this does 
not seriously affect our conclusions. In the external region, then, ordinary 
electromagnetism A^ is unbroken, and the color generator Tg is confined. We 
assume that currents at infinity create a uniform applied Q-magnetic field 
B'^ in the z-direction. We want to know to what degree the fiux is expelled 
from the inner region. 

To study this situation, we need only work in the two-dimensional space 
of gauge symmetry generators spanned by Q and Tg. It is natural to take 
a (the angle by which the unbroken f/(l) is rotated relative to ordinary 



electromagnetism, see ( |2.(j| )) to vary with radius, with a{r) = a^ in the high 
density (inner) region, and a(r) = outside. So in the {Q,Ts) basis. 

To determine what fraction of the magnetic field is expelled from the color 
superconducting region, we must consider what happens at the boundary 
between the two phases. There are three relevant distance scales: 

R the size of the high-density region, 

AR (rl M \-^ ^^^ distance over which a switches from oq to 0, 
the boundary thickness 

A the screening distance for broken/ confined gauge fields 

(2.10) 

We will always assume that i? S> A. We will treat two limiting scenarios: 
"sharp" boundary and "smooth" boundary. The sharp boundary, 5R <^ A, 
corresponds to a sudden step-function-like change in a. This is what we 
would expect to find if there were a first-order phase transition between the 
low and high density regions, with phase boundary thickness less than the 
screening length. The smooth boundary, 5R ^ A, corresponds to a gradual 
change in a. This applies to the situation where there is no first-order phase 
transition between nuclear and quark matter (e.g. for low strange quark mass, 
where we expect no phase transition at all ||2^, |2^, ^), or where there is a 
first-order phase transition with phase boundary thicker than the screening 
length. We will see that the behavior of the magnetic field is quite different 
depending on whether the boundary is sharp or smooth. 

What we are interested in is the behavior of the magnetic field at 
macroscopic distance scales of order R, not at the microscopic scale A. In 
other words, we will always be interested in the unbroken and unconfined 
gauge fields, which obey the Maxwell equations. In order to treat all gauge 
fields in a single formalism it will therefore be convenient to take into account 
screening not by introducing gauge boson masses into the field equations, but 
rather by introducing monopoles and supercurrents to describe the screening 
of the confined and Higgsed gauge fields. By "monopoles" we mean whatever 
gauge field configurations are responsible for terminating magnetic field lines 
of confined gauge fields; we do not require them to be solutions to any classical 
field equations. We will ignore the details of how the supercurrents (and 



"monopoles" ) arrange themselves on distance scales of order A, as all we care 
about is that they screen (terminate) the higgsed (confined) magnetic fields 
incident upon them. 

3 Sharp boundary 



B"^ 



B' 



'± 



r < R 1 r > R 
Higgsed | Confined 

a = ao a = 
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Figure 1: Cross- section of a small region of the sharp boundary, showing the 
location of the X-supercurrents (which ensure that there is no parallel X-fiux 
inside the color superconducting sphere at r < i? — A) and the Tg-monopoles 
(which ensure that there is no perpendicular Tg-fiux outside the sphere at 
r>R + X.) 

In the sharp boundary case, a changes quickly from Oq to at r = i? over 
a boundary region with thickness 6R < A (Fig. |I]). Just inside the boundary, 
in the region R — X < r < R, the X gauge field is Higgsed, so there is a 
density J^ of X-supercurrents that screen out the parallel component of the 
X-fiux, leaving only Q-fiux at r < i? — A. Just outside the boundary, in the 
region R < r < R+ X, the Tg gauge field is confined, so there is a density p^ 
of Tg-monopoles that terminate any perpendicular component of the Tg-fiux, 
leaving only Q-fiux at r > i? + A. 
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To obtain the boundary conditions on the magnetic fields associated with 
the unbroken generators, write the magnetic field in the {Q,Ts) basis ( p.9|) , 

V^^V Uinao^Q + cosao^^ '' ^ ' 



then the Maxwell equations are 

r>R: div5=( Ts), curlS = 0, 
V Pm J 

r<R: div5 = 0, curl5 = ( ~ ^'"""^ "^^ 

cos ao / 



X 



(3.2) 



Since we are only interested in the behavior of the fields on distance 
scales much greater than A, we integrate the Maxwell equations (|3.2|) over 
R — X<r<R + X, and obtain boundary conditions that relate the fields 
at i? — A to those at i? + A. We follow the standard derivation (Ref. pO| , 
sect. 1.5). We find a discontinuity in the normal compoment B^^ due to the 
surface density of Tg-monopoles, and in the parallel component Sif due to 
surface X-supercurrents, 



(3.3) 



B^AR^^){l)-BiiR-x){z::)-^^{i)^ 

II \0/ \smao J \ cosaoj 

From this we immediately obtain the boundary conditions on the fiux: 

5?(i?-A) = -^—B'^{R+X), 

cosao (3.4) 

Bf{R-X) = cosao B^{R+X). 

Thus, just inside the interface we find a B^ whose component parallel 
(perpendicular) to the interface is weakened (strengthened) relative to that 
of B'^ just outside the interface. 

3.1 Solution for the sharp boundary 

Since this is a magnetostatic problem we can write it in terms of a magnetic 
scalar potential $. The potential is associated with the unbroken Q flux 



inside the sphere, and the unbroken Q-flux outside 

B^ = — V$ outside sphere, 
B'^ = — V$ inside sphere. 

Maxwell's equations become 



(3.5) 
V^$ = 0, (3.6) 



with boundary conditions 

$(r — > oo) = B9rcos9, 



or cosao c*^ l-J-' J 

— (i?-A) =cosao^{R+X). 



Expanding in Legendre polynomials (see |Q sect. 5.12) we find the solution 



$ = -B^r cos 9-B^?- cos 6 "' ^"^ {r > R), 

" r2 2 + cos2ao ^ ^' ^ ^ 

(3.8) 

^ ^n ^ 3 COS an / ^N 

$ = -B^r cose r^^ r < i? . 

2 + cos2ao 

In Fig. § we show the resultant field configuration for cosao = 0.5. In the 
real world ao is small, so the field is mostly converted into Q fiux by the 
supercurrents and monopoles, and penetrates the interior. Only a weak field 
is excluded. 

We can check that the solution (|3.8|) makes sense in two limits. (1) cto = 0, 
so Q = Q and Tg = X. The unbroken U{1) in the color superconducting 
region is exactly the same as conventional electromagnetism: the magnetic 
field does not even notice the color superconductor. (2) ao = 7i/2, so 
Q = Ts and X = Q. The region is a conventional superconductor, 
breaking electromagnetism. The magnetic field is completely expelled from 
the superconducting region. 

How large a field must we apply before we begin to destroy the color 
superconductivity in parts of the sphere? The magnitude of the Q-magnetic 
field just outside the sphere is largest at the equator of the sphere, where it 
is given by 

\B<^{r = R+X,e = n/2)\=B^(l+'—^) . (3.9) 

\ 2 + cos^ ao / 
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Figure 2: Locally unbroken magnetic field inside and outside a sphere 
of high- density quark matter, with sharp boundary. This is ( |3.(j| ), with 
cos ao = 0.5. 



This means that the equator of the sphere sees an X-magnetic field given 
by sinao times ( pIOl) . If the color superconductor is Type I, the color 
superconductivity is destroyed in regions of the sphere near its equator when 
the X-magnetic field exceeds the thermodynamic field H^ ~ 10^® Gauss.0 
This requires 

> l±^n^ . (3.10) 

We conclude that because most of the applied field can penetrate the 
superconductor in the form of Q-flux, while only a small fraction of 
the applied field must be excluded, the applied field at which the color 
superconductivity begins to be destroyed is significantly larger than the 
thermodynamic critical field H^ ~ 10^^ Gauss. 



•^If the material is Type 11, flux tubes begin to penetrate when the X-magnetic field 
exceeds i/^, which is of the same order of magnitude. 
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4 Smooth boundary 

For a smooth boundary, 5R 3> A, which means that a changes slowly relative 
to the screening length A. In this case we assume that there is a continuous 
distribution of monopoles and supercurrents. At a given r these produce flux 
in the f/(l) that is (locally) broken. Since the screening length A is small, 
the net flux at a given r must only be in the unbroken U{1). The Maxwell 
equations therefore take the form 



' ' sma;(rjyj y cosa(r, 

curMS(^°^"^)n= j^-^i^«^-^ 



(4.1) 

{ K 

sva.a{r) J j \^ cosal^r 

where B, pm, J are functions of position. Now 

1- f T^ f cosa(r)\} ,. ^/cosa;(r)\ ^ da f — sma(r) 

div <^ S . } ( U = div5 . ; { + Br— ) { 

[ ysina[r) J j ysina[r) J dr \ cos a[r) 

, ( ^ f cos a(r)\^ , ^ f cos a(r)\ ^ _ / — sina(r) 

curl <B { . ) ( ] } = cmlB . ) ( ] + B x Va { ) { 

[^ \sin a[r) J j ysina[r) J y cosa(rj 

'(4.2) 

Substituting these into ( |4.1] ) and taking the (cos a, sin a) component, we find 

that the locally unbroken magnetic field B obeys the sourceless Maxwell 

equations: 

div5 = 0, curl5 = 0. (4.3) 

We conclude that in this case the magnetic field always rotates to be locally 
unbroken, but is otherwise unaffected: it obeys the free Maxwell equations 
everywhere. There is no expulsion at all. (In the smooth boundary case, then, 
the effective critical magnetic field is infinite: an arbitrarily strong magnetic 
field is always rotated to be in the locally unbroken f/(l), and never destroys 
the color superconductivity.) 

One way of understanding this result is to note that the smallness of 
the screening length A means that broken gauge fields are always zero. 
Consequently, when we move from radius r to r + dr, the gauge field is 
immediately projected onto the locally unbroken generator Q{r + dr), which 
differs by an infinitesimal angle from Q{r). As we vary r, we therefore have a 
sequence of such projections. In the limit of an infinite number of projections. 
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each infinitesimally different from tlie last, tlie gauge field is simply rotated 
to follow the locally unbroken generator, 



N 



j™n<=-(ii)=i- (4-4) 



AT^ooJ-l ^ ]\f 



n=l 



This is analogous to the behavior of a sequence of A^ polarizers each at an 
angle ao/N to the previous one. In the limit where A^ -^ oo, an incoming 
photon aligned with the first polarizer exits from the last polarizer with its 
polarization rotated by ao, and with no loss of intensity. 

Finally, it is interesting to ask why the 5R -^ limit of our result for 
the smooth boundary is not the same as our result for the sharp boundary. 
The reason is that in the smooth case we assume A -C SR, so the monopoles 
and supercurrents are always in the region where a is changing. In the sharp 
case, by contrast, we assume 5R <^ A, and there is in effect no such region. 
Taking the 5R —^ limit of the smooth case (A -C 5R) would result in a pile- 
up of currents and monopoles in a shrinking region within which a changes 
rapidly from ao to 0. In this limit, ( [4.3|) is maintained and no flux is excluded 
from the sphere of color superconductor.^ The sharp boundary studied in 
Sect, ^represents different physics, in which A is constant as 5R — > 0. In this 
limit all the monopoles are in the a = region on the confined side of the 
boundary while all the supercurrents are in the a = ao region on the Higgsed 
side and we obtain the boundary condition ( p.4|) and partial flux exclusion 
as in (13T 



5 Consequences for Neutron Stars and Outlook 

If a neutron star features a core made of color superconducting quark matter, 
we have learned that this core exhibits (almost) no Meissner effect in response 
to an applied magnetic field. Even though the Cooper pairs of quarks have 
nonzero electric charge, there is an unbroken gauge symmetry U{1)q, and 

the color superconducting region can support a Q-magnetic field. If the 
boundary layer of the color superconducting region is thick compared to the 

■'Formally, taking the 5R -^ limit while maintaining A ^ 5R yields both monopoles 
and supercurrents concentrated at a point where a — af)/2. If one modifies the right hand 



side of (3.3) to reflect this, the boundary conditions (3.4) becomes B'^{R—\) = B'^{R+X) 



and no flux is excluded. 
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screening length, there is no Meissner effect: the Q-magnetic field within has 
the same strength as the applied Q-magnetic field. (This smooth boundary 
case applies if the color superconducting phase and the baryonic phase are 
not separated by a phase transition; it may also apply in the presence of 
a first order phase transition, if the phase boundary is thick enough.) If 
the thickness of the boundary of the superconducting region is less than the 
penetration length, there is a partial Meissner effect: the Q-Ta.a.gn.eiic field 
within is somewhat reduced relative to the applied Q-magnetic field. Because 
ao is small in nature, the overlap between Q-photons and Q-photons is large, 
and this partial Meissner effect occurs only at the few percent level. 

The physics of magnetic fields in neutron stars with color superconducting 
cores is qualitatively different from that in conventional neutron stars. In 
conventional neutron stars, proton-proton pairing breaks U{1)q and there 
is no unbroken U{1)q. This means that magnetic fields thread the cores of 
conventional neutron stars in quantized flux tubes, within each of which there 
is a nonsuper conducting region. In contrast, the Q- magnetic field within a 
color superconducting neutron star core is not confined to flux tubes. This 
means that, as we discuss below, the enormous Q-electrical conductivity 
of the matter ensures that the Q-magnetic field is constant in time. In 
ordinary neutron stars the Q-magnetic flux tubes can be dragged about by 
the outward motion of the rotational vortices as the neutron star spins down 
3T| , p2| , p!0| , |33| , p4| , and can also be pushed outward if the gap at the proton 



fermi surface increases with depth within the neutron star |H^. One therefore 
expects the magnetic field of an isolated pulsar to decay over billions of years 
as it spins down |^, 0, ^, |3 or perhaps more quickly |^. There is no 



observational evidence for the decay of the magnetic field of an isolated pular 
over periods of billions of years |jTO|; this is consistent with the hypothesis 
that they contain color superconducting cores which serve as "anchors" for 
magnetic field, because they support a Q-magnetic field which does not decay. 
We now estimate the decay time of the Q-magnetic field for neutron 
stars with color superconducting cores, doing the calculation separately for 
cores which are in the 2SC and the CFL phase. To this point, the only 
difference between the two color superconducting phases in their response to 
applied magnetic fields has been a difference of a factor of two in the value 
of ao; because ao <C 1 in both phases, this difference is of no qualitative 
consequence. However, the 2SC and the CFL phase differ qualitatively 
in their symmetries and their low-energy excitations, and can therefore be 
expected to have quite distinct transport properties. As we will see below, 
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the Q-electrical conductivity in the CFL phase is much larger than the 
conductivity in the CFL phase. However, even the "smaller" conductivity of 
the 2SC phase is so large that the timescale for the decay of a Q-magnetic 
field within a color superconducting neutron star core is long compared to 
the age of the universe. 

The characteristic magnetic field decay time due to ohmic dissipation is 



m 



tdecay ~ 4ai?V7r , (5.1) 

where R is the radius of the color superconducting core and a is the Q-electric 
conductivity. (We set h = c = ks = I throughout.) We begin by estimating 
a, and hence the decay time, for a core which is in the 2SC phase. At keV 
temperatures, the dominant carriers are the relativistic electrons and those 
quarks of the 2SC phase which are ungapped, or which acquire gaps so small 
as to be ~ T or less. In the 2SC phase the up and down quarks of one color 
acquire gaps smaller than or of order keV |^ . The strange quarks of all three 
colors — which do not participate in the dominant pairing characterizing the 
2SC phase — can be expected to have gaps which are similar in magnitude or 
even smaller |^.0 To obtain a lower bound on a, we assume that all strange 
quarks, and up and down quarks of one color, have gaps <ti T. These five 
quarks have Q-charges — |,— |,0,l,Oin units of the Q-charge of the electron 



egj \/ g^ + e^/12, which we henceforth take to be just e since ej g <C 1. The 
sum of the squares of the Q-charge of the ungapped quarks is therefore 3e^/2. 



Following [0, the electrical conductivity is given by 



, 2„2„ I ^ , 2„2 



-/i^e Te + - —^11 e Tq , (5.2) 



where Tf. and r^ are the momentum relaxation times for electrons and quarks 
in the plasma, defined and calculated in Ref. |^. For both the electrons 



and quarks, the dominant scattering process contributing to momentum 
relaxation is scattering off quarks. The five gapless quarks of interest yield 



m 



40 T^/^ 

r-^^— l.Sla^^^ (5.3) 



^ The contribution from the quark quasiparticles with gaps A ~ 20 — lOOMeV ^ T 
can be neglected. In the CFL phase, all quark quasiparticles have gaps which are ^ T. 
We discuss this situation below. 
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where Ug = g'^/4:7i and {rnFj^Y ~ 3(7^/i^/(27r^) is the Debye screening mass 
for the gluons, neglecting Mg relative to /i, and where we have assumed that 
T <^ ?^f)-Q Similarly, 

4 t5/3 

r-i^— 1.81a^-— ^, (5.4) 

where ctg = 1/137, and {rrf^Y — 5e^/i^/(127r^) is the Debye screening mass 
for the Q-photon. We have neglected M^, used the fact that the average 
squared Q-charge of the nine quarks participating in screening is 5e^/18, and 
have assumed T <^ Trfj^. Taking a^ ~ 1, we find first that the electrons 
dominate the conductivity, by a factor of about 20, and second that the 
time-scale for the decay of the magnetic field in a color superconducting 2SC 
core of radius R is 

taeca, ~ 3 X 10 years [j^j (^5^1^) (^^l^j [y^) ■ 

(5.5) 
Thus, the magnetic field in the core of a neutron star which is made of matter 
in the 2SC phase decays only on a time-scale which exceeds the age of the 
universe. 

We now turn to a color superconducting core in the CFL phase. In 
contrast to the 2SC phase, the condensation in the CFL phase produces 
gaps ^ T for quarks of all three colors and all three fiavors.Q Thus there 
are no low-energy quark quasiparticle excitations. Furthermore, the CFL 
condensate gives a mass to all eight gluons. In this quark matter phase, 
the degrees of freedom which are most easily excited are neither quarks 
nor gluons. Because of the spontaneous breakdown of chiral symmetry in 
the CFL phase, there are charged Nambu-Goldstone bosons, which would 
be massless pseudoscalar excitations if the quarks were massless. Once the 
nonzero quark masses are taken into account, one finds pseudoscalar masses 
which are small |l^ (in the sense that they are -C A) but which are still large 
compared to T. There is one remaining scalar Nambu-Goldstone excitation. 



^Note that we have worked directly from equations (28) and (39) of Ref. |37), and have 
not used the numerical factor in equation (40) of Ref. |^^, which contains an error [ p8[ . 

^It may be possible for some quark quasiparticles to be gapless even in the CFL phase, 
at densities just above those where the 2SC phase is favored pTJ. In this nongeneric 
circumstance, the CFL phase conductivity would be similar to that of the 2SC phase. 
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associated with the superfluidity of the CFL phase, but this excitation is Q- 
neutral. We thus discover that all possible hadronic excitations which have 
nonzero Q-charge do acquire a gap: their populations in thermal equilibrium 
are suppressed exponentially by factors of the form exp(— A/T), where A is 
either a Fermi surface gap or a suitable bosonic mass. This first of all means 
that the only charge carriers which could contribute to a are the electrons. 
Second, the scattering of these electrons off the positively charged hadronic 
system in which they are immersed vanishes exponentially for T — >^ 0, 
because at low temperatures there are no hadronic excitations off which to 
scatter.]] The electrons can only scatter off other electrons. However, such 
collisions do not alter the total electric current |3^. We therefore conclude 
that although matter in the CFL phase is not a Q-superconductor (it does 
not exclude Q-magnetic field) it is a near-perfect conductor: the resistivity 
drops exponentially to zero as T — > 0. At typical neutron star temperatures 
T ^ A, the density of hadronic excitations, and therefore the resistivity, is 
exponentially suppressed relative to the 2SC phase. As a consequence, the 
decay time for a Q-v[va,gn.eiic field in a CFL core is exponentially larger than 
that (|5.5|) for a core in the 2SC phase, which was already longer than the age 
of the universe. 

It is clear that the Q-magnetic field is rigidly locked in the color 
superconducting core. It cannot decay with time, even if rotational vortices 
move through the core as the spin rate of the pulsar changes with time. 
Rotational vortices do exist in the CFL phase, because the CFL condensate 
spontaneously breaks a global U{1). Instead of dragging magnetic flux tubes 
with them as they move, as occurs in a conventional neutron star, the 
rotational vortices can move freely through the CFL phase because there 
are no Q-flux tubes, only a Q-magnetic field. Thus, as the spin period of 
the neutron star changes and the rotational vortices move accordingly, there 
is no change at all in the strength of the Q-niagnetic field in the core. The 
conclusion is the same in the 2SC phase, but the argument is even simpler 
because in this phase there is no spontaneously broken global f/(l), and 
therefore no rotational vortices. 



^We are describing the conductivity in the hnear response regime. In particular, we are 
assuming that the current is smaU enough that the momenta acquired by the electrons due 
to the current is small compared to the excitation energy for all charged hadronic modes. 
If the current were increased beyond the linear regime, the conductivity would decrease 
and eventually the CFL condensate itself would be destroyed. The nonlinear regime is not 
relevant for our purposes. 
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As we have noted above, the data on isolated pulsars show no evidence 
for any decay of the observed magnetic field even as they spin-down over time 



| I0| . This is consistent with the possibility of a color superconducting core 
within which the field does not decay. However, we must also ask whether and 
how the fact that the observed magnetic fields of accreting pulsars do change 



as they spin up |]T0| is consistent with the possibility of quark matter cores. 
There are several ways in which the surface magnetic field could decrease 
as a neutron star accretes and spins up, even though the magnetic field in 
the core remains constant. One possibility is that the accreting matter may 
bury the magnetic field ^, |T0|. Another possibility is that as the magnetic 



flux tubes in the mantle and crust of the neutron star are pulled around 
by the rotational vortices, the north and south magnetic poles on the star's 
surface may be pushed toward one another, reducing the observed dipole 
field even though the field deep within, in the color superconducting core, 
remains undiminished |^, p4[ . 

The analysis of Section ^ was idealized in three ways relative to that 
appropriate for a neutron star, if there is a first order phase transition 
between baryonic and quark matter. First, we assumed that ordinary 
electromagnetism was unbroken outside the core. This is false: proton-proton 
superconductivity results in the restriction of ordinary Q-magnetic field to 
flux tubes. This means that our derivation of the boundary conditions in 
Section |^ only applies upon averaging over an area of the boundary that is 
sufficiently large that many Q-finx tubes are incident on it from the outside. 
A more microscopic description of the field configuration near the boundary 
would in fact require further work, but this is not necessary for our purposes. 

Second, our assumption in Section |^ of a spherical boundary is 
oversimplified. If there is a first order phase transition between baryonic and 
quark matter, because there are two distinct chemical potentials /i and fie 
there will be a mixed phase region, with many boundaries separating regions 
of quark matter and baryonic matter which have complex shapes [^ . This 
complication of the geometry of the boundary evidently makes a complete 
calculation more difficult than the one we have done in Section ^, although 
if the boundary is thick then the conclusions of Section ^ are unaffected. 
Regardless, the qualitative conclusion that only a very small fraction of the 
flux will be excluded because ao is so small will not be affected by these 
complications. 

Third, the conflguration of Figure 2 cannot, in fact, be attained in a 
neutron star even though it is favored in the sharp boundary case. The 
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core of a newborn neutron star is threaded with ordinary Q-magnetic field. 
Because of the enormous conductivity, the time it would take to accomplish 
the partial exclusion of fiux seen in Figure 2 is exceedingly long ||3^. This 
means that, instead, although the field within the core will be largely Q- 
magnetic field, there will in addition be a small fraction of X-magnetic fiux 
confined in quantized fiux tubes. The sum of the Q- and X-fiuxes adds up to 
the original Q-fiux. Over time, the motion of rotational vortices may move 
the X-fiux tubes around. The much larger Q-fiux, which is not constrained 
in fiux tubes, is frozen as described above. 

We conclude that relaxing the simplifying assumptions which we have 
made would not change our qualitative conclusions. If neutron stars 
contain quark matter cores, those cores will exclude at most a very small 
fraction of any applied magnetic field. Instead, the fiux penetrates (almost) 
undiminished. The only change in the fiux within the color superconducting 
core is that it is a Q-magnetic field, associated with that linear combination 
of the photon and the Tg gluon which is unbroken by the quark pair 
condensate. Most important for neutron star phenomenology, and in 
qualitative distinction from the results for conventional neutron stars, is the 
conclusion that this Q-fiux does not form quantized fiux tubes and is frozen 
over timescales long compared to the age of the universe. 

We are confident that we have not said the last word on the effects of 
color superconducting quark matter cores within neutron stars on magnetic 
field evolution. More work is required to better understand pulsars which 
are accreting and spinning up. We can already conclude that if observational 
evidence were to emerge that an isolated pulsar loses its magnetic field as it 
spins down (in such a way that the field would vanish in the limit in which 
the spin vanishes), this would allow one to infer that such a pulsar does not 
have a quark matter core. 
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